A litmus test for A 
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The critical issue in cosmology today lies in determining if the cosmological constant is the under- 
lying ingredient of dark energy. Our profound lack of understanding of the physics of dark energy 
places severe constrains on our ability to say anything about its possible dynamical nature. Quoted 
errors on the equation of state, w{z), are so heavily dependent on necessarily over-simplified param- 
eterisations they are at risk of being rendered meaningless. Moreover, the existence of degeneracies 
between the reconstructed w{z) and the matter and curvature densities weakens any conclusions 
still further. We propose consistency tests for the cosmological constant which provide a direct ob- 
servational signal if A is wrong, regardless of the densities of matter and curvature. As an example 
of its utility, our flat case test can warn of a small transition from w{z) = — 1 of 20% from SNAP 
quality data at 4-(t, even when direct reconstruction techniques see virtually no evidence for devia- 
tion froin A. It is shown to successfully rule out a wide range of non-A dark energy models with no 
reliance on knowledge of Qrn using SNAP-quaility data and a large range for using 10^ supernovae 
as forecasted for LSST. 



Introduction The discovery that the expansion rate of 
the universe is apparently speeding up from the type-la 
supernovae is arguably one of the most significant events 
in the history of modern cosmology [l|] . Although seem- 
ingly consistent with our current concordance model in 
which the source of the cosmic acceleration, coined 'dark 
energy', takes on the form of Einstein's cosmological con- 
stant, the precision of current data is not sufficient to 
rule out the possibility of an evolving component. A 
huge question mark still exists: whether the canonical 
flat ACDM model is correct or not. If not then we are 
perhaps looking for some dynamical field with a repulsive 
gravitation force, or maybe the field equations are wrong 
on Hubble scales. More bizarrely, this could instead be 
indicating that the Copernican principle is wrong, and 
that radial inhomogeneity is responsible for our confu- 
sion [2j]. Within the FLRW paradigm, all possibilities 
can be characterised, as far as the background dynam- 
ics are concerned, by the dark energy equation of state 
w{z). Unfortunately, from a theoretical perspective we 
have virtually nothing to go on, implying that w{z) could 
really be pretty much anything. Our priority in cosmol- 
ogy today must therefore lie in searching for evidence for 
w{z) ^-\. 

The observational challenge to such ambiguity lies in 
trying to find a straightforward, yet meaningful and suf- 
ficiently general way to treat w^z). This is usually in 
terms of a simple parameterisation of w{^z\, but any such 
functional forms for w(z) are highly problematic because 
they have no basis in a sound physical theory and to 
be flexible require a set of parameters that becomes far 
too large to be meaningful. Furthermore, they are only 
sensitive to deviations from w = — 1 in a highly limited 
function space [3[. To try to reduce the huge arbitrari- 
ness, the space of allowed w{z) models is often reduced 
on physical grounds to w > —1; but if w is only effec- 



tive, parameterising a modified gravity model, say, then 
this might be too limiting. Perhaps the dark energy is 
getting up to all-sorts: until we have a way of predict- 
ing this, prematurely deciding that w{z) must be slowly 
changing or undergoing a sharp transition say, leaves vast 
unchartered blind spots in the possible function space of 
dark energy [3| . Ideally, then, we need a non-parametric 
way of determining the functional form of w(z). 

An alternative procedure is to reconstruct w{z) di- 
rectly from the observables without any dependence on 
a parameterisation of w{z) or understanding of dark en- 
ergy y, la, y ■ Such direct reconstruction methods rely on 
estimating the first and second derivatives of luminosity- 
distance data. Errors on dL{z) alone do not translate 
simply to errors on wiz): dL{z) can deviate from ACDM 
< 1% even though w{z) is fluctuating wildly. The re- 
verse of this was shown by means of striking plots in [7| . 
Furthermore, the precision to which we can estimate the 
matter and curvature densities have acute ramifications 
for estimates oi w{z) |7l.l8l.l9l.lld|. 

Because the quality of the current data is not yet of 
sufficient quality to determine the exact form of w{z) or 
to distinguish between different models with much con- 
fidence, it is difficult to be satisfied that w — —I. Here, 
we demonstrate a method to put the concordance model 
to the test by introducing a simple null test for ACDM 
directly from distance-redshift data. We construct a rela- 
tion ^{z) which has the property of being exactly equal 
to over any redshift range of the reconstruction only 
for w = —I. Critically, it is independent of the matter 
density flm- Using the Fade ansatz introdued in |_53 to fit 
directly to the luminosity-distance curve, we show that 
^{z) can be used to rule out w = — 1 for a wide range of 
types of evolving w{z) for 10^ SN as predicted by LSST 
[ill, with certain test cases oi'w{z) being ruled out using 
only SNAP-like data. We then extend our test to include 



curvature, thereby providing a genuine litmus test for the cosmological constant. 

I 

Reconstructing Dark Energy The dark energy equation of state, w{z), is typically reconstructed using distance 
measurements as a function of redshift. The luminosity distance may be written as 
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which is formally valid for all curvatures, where H{z) is given by the Friedmann equation 

l + w{z') 



H{zY = H^\ a„.(l + zY + Uk{l + zY + ^DE exp 



dz' 



(1) 
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and ^DE = 1 — ^m — ^k- The usual procedure is to postulate a several parameter form for w{z) and calculate dL{z) 
using Eq. ([T]). An alternative method is to reconstruct w{z) by directly reconstructing the luminosity-distance curve. 
It has been shown in [9|? ] dL{z) may be inverted to yield w{z). Writing D{z) — {Hq/c){1 + z)^^di(z), we have 



w{z) 



2(1 + z)(l + nkD^)D" - [{1 + z)2j7feD'2 + 2(1 + z)nkDD' - 3(1 + flkD^)] D' 

3 {(1 + zf [Vik + (1 + z)nj\ D'^ - (1 + nkD'^)} D' ■ 



(3) 



Thus, given a distance-redshift curve D{z), we can reconstruct the dark energy equation of state [? ]. Typically, one 
chooses a parameterised ansatz for D{z), such as the Fade ansatz used in [5|], 



Dl{z) 



jl + z) -a^{l + z) - l+Q 



l3{l + z)+jy/m+2 
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(4) 



one then fits it to the data, and then calculates w{z) from Eq ([3]). Such a reconstruction method is more generic 
than parameterising w(z) directly because we are fitting directly to the observable, and so can spot small variations 
in dL which can translate to dramatic variations in 'w{z). Unfortunately, this also leads to larger reported errors in 
w(z) than calculated by parameterising it directly. It is not clear, however, which errors on our understanding of 
dark energy should be taken most seriously. What is very nice about this method is that, if done in small redshift 
bins, measurements of iy(z) in a given bin is independent of bins at lower z; this is not the case for parameterised w 
methods as they integrate over redshift when calculating dL{z). Both methods suffer from strong degeneracies with 
fife and r2„j, however - see, for example [8|, |9|, |lO| . This vital problem we circumvent using the tests we present below. 
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Consistency Test for ACDM For flat ACDM mod- 
els the slope of the distance data curve must satisfy 
D'{z) = l/-\/17m(l -t- zY + (1 — ^m)- Rearranging for 
flm we have 

l-D'jzY ... 

"'' [il + zY-l]D'{zY' ^' 

Within the flat ACDM paradigm, if we measure D'{z) at 
some z and calculate the rhs of this equation, we should 
obtain the same answer independently of the redshift of 
measurement. Differentiating Eq. ([5]) we then find that 

^(z) = (D"{z) + 3(1 + zYD'{z)[1 - D\zY] 

= for aU flat ACDM models. (6) 

We have used the shorthand ( = 2[(1 -|- zY — 1]- Note 
that this is completely independent of the value of r2„i- 
We may use this as a test for A as follows: take a param- 
eterised form for D{z) and fit to the data. If -Sf (z) = 



lies outside the n-cr error bars, we have n-cr evidence for 
deviations from A, assuming the curvature is zero (see 
below for the generalised test without this restriction). 
If, on the other hand, ^(z) == is a sufficiently good 
fit for all suitable parameterisation we can think of, then 
that is good evidence for A. But the important thing 
is that only one choice of parameterisation has to imply 
^(z) ^ to provide evidence that ACDM is wrong. Ev- 
ery parameterisation has many blind spots - variations 
in w{z) that could never be picked up - so by searching 
through various ensures that the blind spots are illumi- 
nated. 

Testing the test To illustrate the test, and to demon- 
strate that it will work, we fit the Fade ansatz given 
in Eq. Q to the luminosity-distance data for a set of 
simulated data. The derivatives D" and D' are then 
calculated and substituted into Eq. ^. A deviation of 
^(z) = within the error bars of the reconstruction sig- 
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FIG. 1: Consistency test for A. We show Jif{z) for various dark energy models that we could expect to reconstruct from SNAP 
(middle row) and LSST )bottom row) data, which is independent of Q,m- Also shown is the reconstructed iu(z) using the 
Pade ansatz, Eq. Q, together with 1-a errors; this reconstruction assumes the correct choice of flm, and it considerably worse 
without this. For SNAP data, we see that .if (2) always provides evidence for evolving dark energy, except in the case when it 
is evolving slowly, as in the middle column. Note that the Pade ansatz is especially well adapted for an evolving model such 
as that in the fourth column, but is extremely poor for constant w, or where it undergoes a sharp transition (fifth column). 
Indeed, in the transition model, it provides little evidence for w 7^ —1 while .if(z) signals something is happening at A-a. In 
the third row, we see that ^{z) will be able to signal evolving dark energy to high confidence using LSST data. (Note the 
different redshift ranges used for SNAP and LSST.) 



nals the detection of a non-A paradigm. In order to test 
the consistency equation, we select a sub-class of parame- 
terizations encompassing a wide class of of DE models. In 
addition to it; = —1, we consider quintessence- type mod- 
els with w ^constant, slowly evolving w models, and a 
model which undergoes a rapid transition. 

To begin with we use simulated supernovae data which 
we expect to have from the SNAP mission [lj|. The 
results for ^{z) (in red) for w{z) cases discussed above 
from 1000 realizations of SNAP-quality data with with 
2013 SN distributed throughout z — — 1.7 as done in 
[lOJ I with Gaussian noise are shown in the top two rows 

of Fig. m 

We find that with SNAP data we will be able to rule 
out ACDM to more than la in certain redshift ranges for 
most dynamical models, whereas ^{z) = is included 
in the error bars for w — —1 for all z as we would expect. 
The discrepancies in the effectiveness of -Sf (z) for dif- 
ferent models is readily understood. The example w{z) 
in the third and fourth columns represent dark energy 
models in which the equation of state evolves away from 
w = —1 with increasing z. We know that dark energy 
with w < becomes less dominant at earlier times which 
is why its impact on the observables registers to a lessen- 
ing degree with increased z. The reconstructed w{z) for 
each case is shown in green in the top row. Although the 



reconstruction oi w{z) itself appears in some cases to be 
more effective in ruling out A, the results are heavily re- 
liant of the correct choice of flm- Perturbing fl^ used in 
the reconstruction by only 7% from the true value of the 
simulated data renders an incorrect reconstruction and 
includes w = —1 with the la error bars. =S?(z) however 
is unaffected by Qm- Furthermore, we can see that in 
some cases the two methods can be complimentary: for 
the case w — —0.8, we see that the standard reconstruc- 
tion rules out A to >l-a for z < 0.75, while ^{z) rules 
it out for 0.5 < z < 1.3. 

The information that can be derived from a SNAP- 
like data set is not of a sufficient accuracy to rule out 
w = —1 for all of our test cases. We now consider the 
requirements of a future-generation survey to determine 
if this is a robust test of A. We consider a future data 
set of 10^ SN of SNAP-like quality distributed uniformly 
throughout z = — 1 in line with predictions from LSST 
ll|. Figure [T] shows ^{z) for the same set of models on- 
sidered above, in the bottom row. It is clear that w — —1 
is ruled out to at least 2t7 in all cases. Although 10^ SN 
may be regarded as a somewhat optimistic estimate, it 
shows the potential of our consistency relation to test 
for A without knowing f2m at all. Even though the re- 
construction of w{z) will improve with such a large data 
set, its innate dependence on the matter density Qm and 




FIG. 2; Reconstructed w{z) for w = —1 when correct 
fim = 0.275 is assumed (dark green) and when the incor- 
rect value r2„i — 0.2 is assumed (Ught green) with la error 
bars. This is for a data set of 10^ SN as predicted from LSST. 
This demonstrates the short-comings of reconstruction meth- 
ods which rely on precise measurements of the matter density. 
Our consistency test is immune to such pitfalls. 



the method of reconstruction means that the error bars 
although smaller are not the whole story. This point is 
illustrated in Fig. [51 where we contrast the reconstructed 
w{z) for 10^ SN when we assume the correct ilm = 0.275 
and the wrong f2„j — 0.2. 

Including curvature Our consistency test may easily 
be extended to include curvature as follows: with flk ^ 0, 
^{z) becomes: 

^(z) = {(D^D" - [z2(3 + z)(l + z)D'^ + CDD' 

- 3(1 + zfD^]D'} ^k + CD" + 3(1 + zfD'il - D'^) 
= for ACDM. (7) 

In order to formulate a test which is independent of cur- 
vature, we can rearrange this for fi^ and differentiate it 
to yield a (large!) expression whose numerator must be 
equal to zero. Unfortunately, as it involves third deriva- 
tives of D(z), it doesn't appear to be very useful: an anal- 
ysis similar to the above reveals that even 10^ SN is not 
enough. Rearranging for 17^ and comparing at different 
redshifts will also be a test for A. On the other hand, cur- 
vature may be determined directly from D{z) and H{z) 



where h{z) — H{z)/Hq. Given H{z) data from BAO 
measurements or from relative age measurements of pas- 
sively evolving galaxies, we can expect this more general 
test to be just as useful as the simplified version we have 
analysed here by the time SNAP or LSST data is avail- 
able. Note that it involves only second derivatives of 
distance data and no derivatives of the Hubble rate. 



Conclusions We have proposed a simple and direct 
litmus test for the canonical ACDM paradigm. Our flat 
consistency test is shown to rule out w = —1 for a broad 
range of dynamical models of dark energy to at least 
2(7 using 10^ SN as forecasted for LSST while being en- 
tirely independent of the matter density. We have also 
shown that rapid variations in w{z) can be constrained 
extremely tightly by data we can expect very soon. We 
have also proposed a litmus test for the cosmological con- 
stant which is independent of all cosmological parame- 
ters. This may play an important role in determining 
the true nature of the dark energy. 



Acknowledgements We would like to thank Pedro Fer- 
reira, Roy Maartens, Joe Zuntz and Mark Sullivan for 
useful comments. CC acknowledges support from the 
NRF (South Africa). CZ is supported by a Domus A 
scholarship awarded by Merton College. 



data using ^k = {[H{z)D'{z)Y - l}/[HQD{z)f [9|. Sub- 
stituting this into Eq. JT]) gives an observational litmus 
test for A which is independent of all cosmological density 
parameters: 

^ge„(z) - {CD^D" - [z^i + z){l + z)D'^ 

+ CDD' - 3(1 + zfD^]D'] h{zf + CD 
-^•i{\ + z)[?>{l + z)D^-z^{z + i)]D\ (8) 
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